Impurity-induced diffusion bias in epitaxial growth 
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We introduce two models for the action of impurities in epitaxial growth. In the first, the in- 
teraction between the diffusing adatoms and the impurities is "barrier" -like and, in the second, it 
is "trap" -like. For the barrier model, we find a symmetry breaking effect that leads to an overall 
down-hill current. As expected, such a current produces Edwards- Wilkinson scaling. For the trap 
model, no symmetry breaking occurs and the scaling behavior appears to be of the conserved-KPZ 
type. 
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The effect of impurities on growth rate and morphology 
is a classic topic of crystal growth theory The most 
thoroughly studied case is the step flow growth of a vic- 
inal surface, when the (immobile) impurities pin the ad- 
vancing steps and thus lead to step bunching |2||] . These 
theories are mesoscopic rather than microscopic in na- 
ture, in the sense that they describe the interaction of 
preexisting steps with discrete impurities |]] or an impu- 
rity concentration field |2| . 

The advent of modern crystal growth techniques aimed 
at manufacturing layers of atomic scale thickness, no- 
tably molecular beam epitaxy (MBE) Q, has lead to a 
renewed appreciation of the fact that small concentra- 
tions of impurity atoms on the growing surface can dras- 
tically influence the growth kinetics. A particularly strik- 
ing aspect of these recent results is that the impurities 
may cither lead to a deterioration of the growth quality 

- as would be expected according to the classic view |IJ] 

- or, conversely, they may play the role of surfactants in 
stabilizing smooth, layer- by-layer growth |^| . 

A clear example of the former type is the effect of hy- 
drogen on the MBE of silicon j|,|7| . It was observed that 
the presence of H in the growth chamber during Si de- 
position on Si(OOl) leads to a decrease of the epitaxial 
height H , at which epitaxy breaks down and the growth 
becomes amorphous, proportional to the logarithm of the 
partial pressure of H . The experiments ruled out the 
hypothesis that the breakdown of epitaxy might be due 
either to an increased coverage of H at the interface ||] 
or to its incorporation into the bulk M. Rather, it was 
concluded that the hydrogen greatly speeds up the 
development of surface roughness due to a reduction of 
the diffusion length of Si adatoms || . 

The modification of the diffusion properties of the 
adatoms appears to be the most significant effect of 
the impurities also when they act as surfactants pO| , pl| , 
though the nature of the modification — for example, 
whether the diffusion length is increased or decreased — 
depends on the chemical species in a complicated way || . 
This sensitivity to atomic details is rather unexpected, 
and calls for the development of models which are more 



microscopic than previous approaches ^JsJ . 

In this Letter, we introduce two models for the ac- 
tion of impurities in epitaxial growth. Rather than at- 
tempting a detailed description of some particular mate- 
rial, our aim is to define a 'minimal' model in which the 
consequences of the impurity-adatom interaction on the 
large scale morphology of the surface can be clearly eluci- 
dated. The study of oversimplified models [ p~2)J]~^ ] of 'ideal 
MBE' |Q has previously been very successful in clarify- 
ing the universality classes for kinetic roughening jl4| [l^] 
and morph olog ical instability |T[|[l7],|l8| in the absence of 
impurities |19|] . 

Our models reproduce the sensitive dependence 
on microscopic details mentioned above: Using two 
equally plausible microscopic interaction mechanisms - 
of 6arner-type and of trap-type, respectively - we find 
that for the barrier model the impurities neutralize the 
destabilizing effect of step edge barriers (ll],[ll| and thus 
lead to smoother growth f20|| , while for the trap model 
the asymptotic morphology remains unaffected. In terms 
of the coarse-grained continuum description of the sur- 
face jl5],[l9| we are able to trace the difference between 
the two models to the fact that the barriers modify the 
symmetry of the surface diffusion process, while the traps 
do not. 

The models proposed in this Letter have three main 
ingredients (see Fig. |l|). 

i) Deposition and Diffusion: For simplicity, we con- 
sider a one-dimensional discrete substrate. Material is 
randomly deposited at a rate F. The deposition occurs 
in a solid-on-solid (SOS) manner, i.e. deposition at a po- 
sition i implies that the surface height h(i) is increased 
by one unit. Every atom which has only one occupied 
neighbor (namely, in the layer below) is considered a mo- 
bile adatom which diffuses at a rate D; atoms with more 
bonds are immobile. When an adatom diffuses to a step 
edge from above we implement an additional energy bar- 
rier pl} | by accepting a diffusion move down the step only 
with probability p = exp(-Es), where the barrier energy 
E$ is measured in units of ksT. In the absence of impu- 
rities such step edge barriers are know to lead to unstable 



1 



growth gum. 
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Fig. 1 

FIG. 1. Schematic representation of the processes present 
in our models. Material is deposited at a rate F and sin- 
gle atoms (gray circles) diffuse at a rate D. At step edges, 
adatom motion down the step is accepted with probability 
p = exp(— Es), where Es is the step edge barrier. Impuri- 
ties (squares) are deposited at a rate F' and are not allowed 
to diffuse. When an atom is deposited over an impurity, the 
latter evaporates and is replaced by the deposited atom. Im- 
purities can interact with the adatoms through two distinct 
mechanism. In the barrier-like interaction, an adatom trying 
to diffuse on top of an impurity will have its move rejected. In 
the trap-like interaction, the same adatom would have been 
trapped, i.e., it would swap positions with the impurity and 
stop diffusing. 

ii) Impurities: Based on the experimental observations 
described previously, we assume that there is a flux F' 
of impurities onto the growing surface. We restrict our 
study to the limit in which the diffusion rate of the im- 
purities is much smaller than the diffusion rate of the 
adatoms, so that the impurities can be considered im- 
mobile. We also assume that impurities evaporate from 
the surface at a rate that keeps 9i, the impurity cover- 
age, approximately constant. This is done by removing 
the impurity whenever a new atom is deposited on top 
of it. An important consequence of this rule is that the 
average lifetime of an impurity at a given site equals the 
monolayer deposition time. 

Hi) Interactions: Concerning the interactions between 
the impurities and the diffusing adatoms, we introduce 
two alternative models: the barrier and the trap model. 
In the barrier model, an adatom trying to diffuse onto 
a site occupied by an impurity will have its attempted 
move rejected. In the trap model, the same adatom will 
move on top of the impurity and then will swap positions. 
The end result is that it will no longer be able to diffuse 



because it has (at least) two "chemical bonds" 
atom below and to the impurity. 

Fig. 2 
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FIG. 2. Plot of the total surface width as a function of time 
(or coverage) for barrier-like impurities. The results shown 
were obtained for L = 1000 and D/F = 10 4 . Avera ges were 
taken over 50 runs. In the absence of step edge barriers or im- 
purities W(t) diverges with an exponent /3 » 1/3. When step 
edge barriers are introduced, the instability in the growth pro- 
cess leads to an effective exponent close to one at intermediate 
times and the random deposition value f3 « 0.5 at long times. 
The presence of even a small amount of impurities leads to a 
significantly smoother surface. The lines are plotted as guides 
to the eye, and have slopes 0.24 (dotted), 0.33 (dashed), and 
0.50 (full). 

In our simulations we focus on the exponent (3 de- 
scribing the increase of the surface width W(t, L) = 
((h — {h)) 2 ) 1 / 2 ~ t& in the early time regime t <C L z ; 
here L denotes the system size and the dynamic expo- 
nent is z = 1/(1 — 2/3) for the class of (one-dimensional) 
models considered in this paper |l||. To put the results 
into perspective, we will compare them to the predic- 
tions of the appropriate continuum equations for the 
coarse-grained height function h(x,t); as usual, the av- 
erage height (h) = Ft will be subtracted. In the present 
context the following equation suffices l|l4|-|iT 



dh 



= vV 2 h 



X 



V 2 (V/i) 2 -n(V 2 ) 2 h + n, 



(1) 



where the stochastic force r)(x, t) models the shot noise 
in the beam, and can be taken to be Gaussian with zero 
mean and covariancc 



(r)(x, t)r)(x',t')) = FS(x - x')5(t - t'). 



(2) 



The first term on the right hand side of ([!]) arises from the 
gradient expansion of an inclination-dependent, growth 
induced surface current |p5|,^8|. When it is present, it 
dominates the large scale morphology: For v > one ob- 
tains kinetic roughening of the Edwards- Wilkinson (EW) 
universality class [011,0 with f3 = 1/4, while for v < 
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Fig. 3(a) 



the growth is unstable and a mound morphology is ex- 
pected to develop ||l7|jl£|^| . If, for reasons of symme- 
try (see below) v — 0, the second, nonlinear term be- 
comes important, and changes the roughening exponent 
to [3 = 1/3 |l5|Ji4| , |l(| (the "conserved Kardar-Parisi- 
Zhang" universality class p4j| ). Finally, in many cases 
the growth-induced coefficients v and A are small, and 
the early time behavior is dominated by the third term 
in (0) , which arises from equilibrium surfa ce diff usion fl25fl 



and leads to a (transient) value (3 = 3/8 [112 13 



We consider first the case of a nonzero step edge bar- 
rier in the absence of impurities; then v < |l5|Jl8| ] and 
one expects asymptotically unstable growth. As shown 
in Fig. ||, the instability sets in after an initial power 
law transient which terminates at about 100 monolayers. 
Later on, the destabilizing effect of the step edge bar- 
rier leads to wavelength selection and mound formation 
with a very rapid growth of the surface width. Once large 
slopes have appeared on the surface, there is hardly any 
transfer of matter between the different mounds and the 
exponent (3 reaches the limiting value (3 — 1/2 character- 
istic of random deposition [^6| . 

When we introduce impurities of the barrier-type a 
striking change occurs. As is visually apparent from Fig. 
||, the presence even of small amounts of barrier impu- 
rities leads to a significative decrease of the interface 
width. The exponent takes the value (3 — 0.24 ± 0.03, 
consistent with EW universality. The natural interpreta- 
tion is that the impurities have caused the coefficient v 
in Eq. ([[]) to change sign, from destabilizing (y < 0) to 
stabilizing (v > 0). The value of v can be directly ascer- 
tained by measuring the average surface current for tilted 
substrates jl8| . The results, shown in Figure |^(a) , con- 
firm our interpretation: Even a small flux of impurities 
(F'/F = 0.05) leads to a sizable positive value of u, both 
in the presence of a step edge barrier and for E$ = 0. 

The reason for the change produced by the barrier im- 
purities can be understood as follows. As is well known 
[|J|Jl8||, the step edge barrier leads to an up-hill current 
because adatoms are rejected when trying to go down 
step edges and become integrated in the bulk when reach- 
ing an ascending step. This difference leads to an average 
current towards the up-step which destabilizes the sur- 
face, as described earlier. To visualize the effect of the 
impurities on this process, let us consider a step train 
moving from left to right; cf. Fig [|(b). An impurity can 
be deposited anywhere on a given terrace, so we can say 
that on average it is deposited in the middle between the 
two steps. However, as more material is deposited, the 
step edge to the left of the impurity advances towards 
it. On the other hand, the step edge to the right of the 
impurity moves away from it. Thus, the distance to the 
step to the left of the impurity is typically smaller than 
the distance to the right. Since the current away from 
the impurity on each side is proportional to the material 
deposited there (and thus to the length of that part of 
the terrace) we see that an average down-hill current is 
generated. 



0.005 



o.ooo 



-0.005 



-0.010 



-0.015 




//' 



Vr 



Fig. 3(b) 



FIG. 3. a) Plot of the current j(m) as a function of the tilt 
m of the interface. The tilt is prescribed numerically through 
helicoidal boundary conditions. Our results make it clear that 
the barrier-like impurities lead to a negative current even in 
the presence of a step edge barrier, while the trap impurities 
leave the sign of the current unaffected. We used a smaller 
value of the step edge barrier because of the numerical diffi- 
culties in calculating the current (see Ref. |l8|). b) Schematic 
representation of the effect of a barrier-like impurity on the 
motion of a step train. On average an impurity will be de- 
posited on the middle of a terrace. However, while the upper 
step moves towards the impurity, the lower step moves away 
from it, so that in fact < l r , where li (l r ) is the length of the 
portion of the step to the left (right) of the impurity. Since 
the up-hill current is proportional to li and the down-hill cur- 
rent is proportional to l r , we will have an average down-hill 
current which leads to a positive v coefficient and a stable 
interface. Note that this mechanism can only be effective if 
during the time to deposit one monolayer the impurity will 
have been removed from its position; otherwise, the impurity 
would pin the step and destabilize the surface. 
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In contrast, the trap impurities do not seem to signifi- 
cantly change the dynamics of the growth process, apart 
from an increase of the prefactor of the width (Fig. |^); 
certainly they are not able to suppress the destabilizing 
effect of the step edge barriers. This is confirmed by a 
measurement of the surface current, which remains up- 
hill in the presence of traps, as shown in Fig. ||(a). 

The traps do not induce a surface current because they 
cannot bias the diffusion of adatoms: The trap impurity 
makes itself felt only when the atom has already jumped 
onto it, and therefore no longer participates in the mass 
transport on the surface (when the impurity disappears 
due to the deposition of an additional atom, the trapped 
adatom remains immobile). 

This point is brought out more clearly by considering 
the trap model without step edge barriers. For E$ = 0, 
the pure model (F' = 0) has a symmetry which forces 
v = in Eq. (|): For any local environment, the proba- 
bility of a mobile adatom to jump to the right is equal to 
that for a jump to the left. Since this is true irrespective 
of the overall surface tilt, no growth-induced current can 
exist [^7| . With v — the behavior of is dominated by 
the second, nonlinear term, and one expects [3 = 1/3. Our 
simulations lead to a value of j3 = 0.33 ± 0.03, indicating 
that this symmetry is preserved by the trap impurities. 

The increase in the prefactor of the width with increas- 
ing trap concentration can be interpreted as a decrease of 
the diffusion length £jy: Indeed, it can be shown ]2£| that 
the prefactor scales as ^~ 4d /( 10+d ) f or a d-dimensional 
surface. From the data shown in Fig. |^ we therefore es- 
timate that an impurity flux F'/F — 0.1 decreases the 
diffusion length by almost a factor of three. 

Fig. 4 
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FIG. 4. Plot of the total surface width as a function of 
time for trap-like impurities, under the same conditions as in 
Fig. |^. In this case, the effect of the impurities is not as strik- 
ing as for the barrier-like impurities. Even in the presence of 
step edge barriers, the main effect seems to be an increase in 
the prefactor of the width. The straight lines are plotted as 
guides to the eye, and have slope 0.33. 

Finally, one may ask why none of our models shows 



an impurity-induced growth instability of the kind con- 
sidered in the classic theories of step bunching We 
believe that this is due to the fact that the lifetime of 
an impurity at a given position is, in our models, fixed 
to be of the order of the monolayer deposition time. In 
terms of the conventional step flow picture |l],|| it is evi- 
dent that impurities can effectively pin steps only if they 
remain at a position much longer than the time required 
for a step to pass over a terrace. Thus, it appears impor- 
tant to consider models with a variable lifetime for the 
impurities. 

We acknowledge stimulating discussions with M. Rost, 
M. Schimschak and P. Smilauer. 
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